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ABSTRACT 

We examine various implications from a dynamical and chemical model of globular 
clusters (GCs), which successfully reproduces the observed abundance patterns and 
the multiple populations of stars in these systems assuming chemical enrichment from 
fast rotating massive stars. Using the model of Decressin et al. (2007) we determine 
the ratio between the observed, present-day mass of globular clusters and their initial 
stellar mass as a function of the stellar initial mass function (IMF). We also compute 
the mass of low-mass stars ejected, and the amount of hydrogen ionising photons 
emitted by the proto globular clusters. Typically, we find that the initial masses of 
GCs must be ^ 8-10 times (or up to 25 times, if second generation stars also escape 
from GCs) larger than the present-day stellar mass. The present-day Galactic GC 
population must then have contributed to approximately 5-8% (10-20%) of the low- 
mass stars in the Galactic halo. We also show that the detection of second generation 
stars in the Galactic halo, recently announced by different groups, provides a new 
constraint on the GC initial mass function (GCIMF). These observations appear to 
rule out a power-law GCIMF, whereas they are compatible with a log-normal one. 
Finally, the high initial masses also imply that GCs must have emitted a large amount 
of ionising photons in the early Universe. Our results reopen the question on the 
initial mass function of GCs, and reinforce earlier conclusions that old GCs could 
have represented a significant contribution to reionise the inter-galactic medium at 
high redshift. 
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1 INTRODUCTION 

Although for a long time thought to be among the most sim- 
ple stellar systems, globular clusters (hereafter GCs) have 
been subject to intense studies both observationally and 
through theory and simulations. These include for exam- 
ple detailed work on the stellar content of GCs and on 
chemical abundances of GC stars, searches for viable proto- 
GCs, studies of dynamical effects on massive star cluster 
evolution, cosmological simulations of their formation, es- 
timates of their contribution to cosmic reionisation, and 
other related Galactic a nd extragalactic astrophysical top- 
ics (see e.g. r eviews by Gratton. Sncdcn & Carrctta '2 0041: 
Brodie fc Strader.,200& : .PiottQ.2009. : .Boilv..20ia : .Elmegreeiil 
2010bl ). 



Despite these studies many open questions remain, 
concerning both GCs as individual objects and as a col- 
lective population. For example, the origin of Galactic 
halo stars and the contribution of GCs to this popula- 



tion, are still unclear (see e.g. iHut fc Diorgovskil 19921: 
IParmenti er fc Gilmord l2007l : iBell et all l2008l : iBolev et all 
2009). Similarly, the shape of the globular cluster ini- 
tial mass function (GCIMF), the nature of the present- 
day globular cluster mass function, and the processes 
and the timescales responsible for tr ansforming the for- 
mer into the latter, a r e debated (see Fall fc Z hang 2001; 
Vesperini fc Zepj|2003l : iParmentier fc Gilmord f2005. 2003; 



Elmegreenll2010al '). Also, first steps are being made in or- 
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der to understand GC format i on in cosmolog ical simu- 
lations llBromm fc Clarkd l2002l: iKravtsov fc Gnc din 20qi; 
iBolev et al.ll2009l : lGriffen et al.ll2010 !). Finallv lRicottil (|200a ') 
has shown that GCs emit enough ionising photons to 
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reionise the Universe, provided their escape fraction /esc 
is of order unity. Examining this question is also of inter- 
est in the present context, where it appears that galax- 
ies found so far in deep surveys are in sufficient to r e ionise 
the inter-galactic medium (see e.g. iBunker et al. I I2OIOI : 
lOuchi et al'l I2OO9I : iMcLure et all I2OI0I : iLabbe et al.l l2010l {. 
and where the main sources responsible of cosmic reion- 
isation, presumably faint, low-mass ga laxies, below the 
current detection limits (cf. ,Choudhurv fc Ferraral l2007l : 



Decressin. Charbonnel k, Meyne^ l2007l: [P'Ercole et all 



IChoudhurv. Ferrara &: Gallerani 20081 '). remain thus to be 
identified. 

A major paradigm-shift has occurred recently in the 
GC community, that sheds new light on these key ques- 
tions. Indeed, detailed abundance studies of their long-lived 
low-mass stars made possible with 8-lOm class telescopes, 
together with high-precision photometry of Galactic GCs 
performed with HST, have revolutionised our picture of 
these stellar systems. It is now clear that individual GCs 
host multiple stellar populations, as shown by their differ- 
ent chemical properties and by multimodal sequences in the 
colour-magnitude diagrams ( Bcdin ct al. 2004; Piot to et al.l 
I2OO7I : iMilone et al. 2008, 2010: Villanova ct al. 20l3) In- 
deed, although nearly all GC^J appear to be fairly homo- 
geneous in heavy elements (i.e., Fe-peak, neutron-capture , 
and alpha-elements, se e e.g. I James et al.|[2004l : ISnedenll2005l : 
ICarretta etlo] bOOQal ). they all exhibit large star-to-star 
abundance variations for light elements from C to Al that 
are the signatures of hydrogen-burning at high tempera- 
ture imp lanted at birth in their long-lived low-mass stars 



(see e.g. Gratton et al. 200ll : Gratton. Sneden fc Carrettal 



20041: ISne don 2005; 'Prantzos. Charbonn el fc lUadisI l20ol [ 
Carretta et al. 2009b, : Charbonnel 201 0^. and reference 
therein). In fact, the so-called O-Na anticorrelation is ubiq- 
uitous in Galactic GCs, and is now accepted as the decisive 
observational criterion distinguishin g bona fide GCs from 
other clusters (|Carretta et al.ir2010lj ). 

The current explanation for these chemical patterns is 
the so-called "self-enrichment" scenario that calls for the 
formation of at least two stellar generations in all GCs 
during their infancy: The first generation stars were born 
with the proto-cluster original composition, which is that 
of contemporary field halo stars, while the second genera- 
tion stars formed from original gas polluted to various de- 
grees by hydrogen-burning processed material ejected by 
more massive, short-lived, first gen eration GC stars. Details 
and re ferences can be found e.g. in lPrantzos fc Charbonn"eil 
l|2006l ). who discuss the pros and cons of two versions of this 
"self-enric hment scenario", invoking either massive AGB 
stars (e.g. |Co ttrell fc Da Costal Il98ll : IVentura et alj I2OOII : 
IVentura fc D'Antona .2009i^. or fast rotatin g massive stars 
(e.g. Decressin. Charbonnel fc Mevnetll2007^ as polluters. 

Whatever the actual polluting stars, an immediate 
consequence of this scenario is that, in order to reproduce 
the present proportion of first to second generation stars 
with acceptable values of the polluters IMF, the initial 
stellar masses of GCs must have been considerably larger 
than their present-day value jPrantzos fc Charbonneill2006l : 



^ With th e notable exception of oj Cen, M22, and M54 (see e.g. 
Da Costa et al. 2009; Johnso n fc PiIachowskill201Cl : IsTeeel et alj 



20081 : iDecressin et all |201G| : icarretta et al.H2010bh . How- 
ever, most of the extra-galactic studies have not yet 
incorporated this revised picture, or have not yet ex- 
plored the resulting implications. Furthermore, the recent 
discovery of stars with signatures characteristic of 2"^^ 
;eneration GC stars among the metal-poor halo population 
Carretta et all l2010bl : iMartell fc Grebell I2OI0I ) sheds new 
light on the amount of low-mass stars ejected from GCs 
and on the initial mass function of these clusters, as we 
shall show below. 

In the present paper, we explore several consequences 
of this ne w paradigm, based on the model that was de- 
veloped by iDecressim Charbonnel fc Mevnetl (|2007l . here- 
after DCM07) to describe the early chemical and dynam- 
ical evolution of GCs. In this model fast-rotating massive 
(A/ > 25 Mq) stars are responsible for the GC pollu- 
tion. The model successfully explains the observed abun- 
dance patterns of present-day GC stars, and has also 
bee n tested with N-body and hydrodynamical simulations 
(see Decr essin. Baumgardt fc Kroupall2008l : IDecressin et al.l 
[2010). Its main assumptions are briefly described and sum- 
marised in fj2l Within this framework, we constrain the re- 
lation between the initial and the present stellar mass of 
GCs (33, as well as the contribution to the stellar halo 
{^^, taking the recent observational ide ntification of sec- 
ond generation stars in the G alactic halo (|Martell fc Grebell 
I2OI0I : ICarretta et al.l [20101)1 ) into account. Implications on 
the GCIMF are derived in SjS) Finally, we derive in ij6] a 
well-defined ionising photon production rate for proto-GCs, 
taking all the detailed observational constraints from nearby 
GCs into account, and estimate their contribution to cosmic 
reionisation. Our main conclusions are summarised in fJT] 



THE ADOPTED CHEMICAL AND 
DYNAMICAL EVOLUTION MODEL 



200' 



'TI: ICarretta et al.ll2010al . and references therein). 



DCM07 and IDecressin et al.1 l|2007l ) have shown that the O- 
Na anticorrelation observed in GC stars can be explained if 
a second generation of low-mass stars form from the ejecta 
of first stellar generation fast rotating massive stars mixed 
with some original interstellar material. In their model, the 
first generation forms the full mass spectrum of stars de- 
scribed by a power-law IMF at the high end and a log- 
normal below 0.8 Mq. The second generation of "polluted" 
stars is assumed to form only low-mass stars, following 
the same log-normal IMF (see below). The model allows 
for dynamical cluster evolution, and more specifically for 
the evaporation of stars due to primordial gas expulsion 
driven by supernovae as well as for long-term dynamical 
processes as described bv 'D ecressin. Baumgardt fc Kroupal 
(2008); Decressin ct al. (2Q1^. 

The main free parameter of the model is the IMF 
slope X above 0.8 MqoI the first stellar generation, the 
low-mass IMF being set for both the first and second 
stellar populatio ns to the pr esent-da y mas s function ob- 
served in GCs (|Paresce fc De March! [2OO0I) . Second gen- 
eration low-mass stars are then formed from the mass 
of slow wind ejecta /sw predicted by the stellar evolu- 
tion models of IDecressin et al] l|2007h and after dilution of 
this material with interstellar gas. The parameter d de- 
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scribing this diluti on is inferred from the observed Li-Na 
anticorrelation (see Decressin. Charbonnel fc Mevnetll2007l : 



ICharbonnel fc DecressinI 1201]} ). We adopt d = 1.15 from 
DCM07 as our standard value and comment on the (rel- 
atively weak) dependence of our results on this parameter. 

Allowing for the escape of a fraction of first generation 
stars, the model then predicts the relative number of 1°* and 
2"'' generation stars, as well as detailed abundances ratios 
of these stars, which successfully reproduce observed abun- 
dance patterns and anticorrelations (see DCM07). 

The fraction of "unpolluted", pristine 1°* generation 
long-lived stars fp still present today in GCs can be deter- 
mined observationally from the di stribution of stars along 
the Q -Na anticorrelation (see e.g. IPrantzos fc Charbonnell 
I2OO6I ). In the nomenclature of DCM07 one has fp = 
where nj^ and are the number of I''* 
and 2""* generation low-mass, long-lived stars, respectively. 
Obser vations of the O - Na ant icorrelation in a large GC sam- 
ple bv lCarretta et all l|2010bl ) provide a median value (±68 
% CL) of fp = O.SS^Q Qg both for the total sample and for 
the lowest metallicity ([Fe/H] < —1) hence oldest subsam- 
ple. The semianalytical model of DCM07 predicts fp as a 
function of the IMF slope, the dilution parameter d, and 
the fraction e]^ of low-mass, 1°* generation stars being lost 
from the cluster due to dynamical processes (see Eqs. 20, 23 
of DCM07). We can therefore invert this problem to deter- 
mine, for each value of the IMF slope, the lost stellar mass 
fraction e]^ from the observed value of fp. With this at 
hand, all the properties of the two stellar generations mixed 
within the GC can be determined (see DCM07). Here we are 
in particular interested in the relation between the present- 
day, observed stellar mass and its total, initial value, as well 
as in the mass of stars ejected from the cluster. These are 
derived below. 

To do so, we generalise the dynamical evolution sce- 
nario discussed in depth by DCM07 that allows for mass 
loss from the cluster (as described by Bll > in their "Sce- 
nario 11") due to mass segregation and evaporation of stars. 
We consider the IMF slope x above 0.8 Mq as a free pa- 
rameter, and we determine the allowed values of x from the 
observed value of fp given above. However, it is understood 
that most observations for proto-GCs indicate an IMF slope 
close to Salpeter (x = 1.35) in this mass range (see e.g . 



Chabrier" 2003 ; 'Pe Marchi. Parosce fc Portegie s Zwartll201oi : 
Bastian. Covov fc Mover 2010). While DCM07 assume that 



all second generation stars are retained within the cluster 
(cll = in their notation), we will subsequently relax this 
assumption, motivated by recent findings of some chemi- 
cally polluted, second generatio n stars in the Galactic halo 
l|Carretta eralll2010bl : IMarteU fc GrebelUMlOiy 

We now follow the semianalytical model of DCM07 and 
their notation. The current mass of 1"* generation long-lived 



(i.e. low-mass) stars, M^, in a GC is M^j^ = fp x Mobs, 
where fp is the fraction of 1°' generation stars determined 
from observations (cf. above), and Mobs is the "observed", 
current stellar mass of the globular cluster, excluding stellar 
remnantfl Allowing for dynamical mass loss of stars from 



^ Following standard stellar evolution, stellar remnants consti- 
tute approximately 30% of the total cluster mass after 12 Gyr. 



the cluster, the total (initial) stellar mass of the cluster can 
be written as 
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where fl2 stands for the fraction of stellar mass forming low- 
mass stars (in the I''* generation), i.e. the mass fraction of 
the IMF found at < 0.8 Mq , and where ej^, is the fraction of 
low-mass stars from the I''' generation having escaped from 
the cluster during its history. The total mass of ejected low- 
mass stars is then 
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where we also allow for a fraction e^L of low-mass stars of 
the second generation to escape (see Q. 

In the dynamical GC scenario discussed by DCM07, e}^ 
can be determined from the observed fraction fp of 1°* gen- 
eration stars for a given slope x of the massive stars IMF, 
assuming a value for the global dilution factor d (see their 
Fig. 40 Since /ll depends only on the IMF, it is straight- 
forward to compute the relation between the observed and 
the total initial mass of GCs (Eq. [T| and the amount of 
stellar mass ejected (Eq. [2]). 



3 THE INITIAL MASS OF GLOBULAR 
CLUSTERS 

The ratio between the initial and present-day stellar mass 
computed in this manner is shown in the left panel of Fig. [1] 
Blue lines show the case of ei^ = as in DCM07, the green 
and magenta curves when accounting for = 0.43 and 
0.65 respectively (see below). As can be seen, the steeper 
the IMF, the higher the ratio m between the initial and cur- 
rent mass. This is the case since in the present framework 
massive stars are responsible for the chemical pollution of 
the GC that leads to the formation of second generation 
low-mass stars. Hence for a steeper IMF, a larger total mass 
of stars is required to compensate the relative decrease of 
massive stellar polluters in order to reproduce fp. In this 
way, the same total mass of massive star ejecta incorporated 
into the second generation stars can be produced. As men- 
tioned above, all results shown here are computed adopting 
d = 1.15 for the dilution parameter. A stronger dilution 
of the ejecta from the rapidly rotating massive stars (i.e. 
increasing d) would imply lower values of Mini /Mobs and 
Mjjc/Mobs, since more material is then available to form 
the 2"^* generation stars. In practice, changing d by a factor 
2 (3) around our "standard" value implies changes of ~ 40- 
50 (100) % in the initial and ejected masses, comparable to 
the differences between the different cases illustrated in Fig. 

m 



This fraction may, however, depend also on th e dynamical evolu- 
tion of the cluster l lKruiissen &: Lamer£ir2008l . cf.). 
^ One obtains = 1 — 
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from their Eqs. 



and 23, for the assumptions of scenario II, but relaxing the 



hypothesis of 
stars. 



: 0, i.e. allowing also for loss of 2" generation 
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Figure 1. Left: Ratio between the initial and present-day mass of globular clusters as a function of the massive star IMF slope (1.35 is 
the Salpeter value) and for /p = O.SS^Q gg (boundaries shown by dotted lines). The blue lines correspond to the case where no second 
generation stars escape from the cluster (e^^ = 0) , while green and magenta lines show the values predicted when allowing also for the 
evaporation of 2nd generation stars with fractions = 0.43 and 0.65 respectively. Right: Mass ratio between the ejected low-mass 
stars (total including 1*^* and 2"^^ generation stars) and the present-day mass as a function of the massive stars IMF slope. As for the 
left figure, blue, green, and magenta lines shown the predictions for = 0, 0.43, and 0.65 respectively, and for the observed range of 
/p. The black lines show the contribution of 2"'^ generation stars for the latter two cases. The shaded area indicates the region where 
the present-day GCs would overpredict the amount of halo stars (assuming 2% of the stellar halo mass in GCs). 



For a Salpeter IMF as assumed in DCM07, the initial 
cluster mass is found to be ^ 8-10 times larger than the 
current (observed) mass when no second generation stars 
are lost, as seen in Fig. [T] Considering massive AGBs as 
potential GC polluters (instead of massive stars), and as- 
suming that none of the seco nd g eneration stars is lost , 
iPrantzos fc Charbonnell l|2006h and ICarretta etall l|2010bl ) 
also found that the original cluster population should have 
been larger than the current one by ~ one order of magni- 
tude for a Salpeter IMF. This agreement is dictated by the 
amount of initial mass of polluters needed to provide enough 
material for the second stellar generation. If we attribute 
the recently observed 2"^^ generation stars in the halo to the 
present population of GCs (cf. 311) we must allow for a loss 
of 2""* generation stars (i.e. > 0), which implies even 
higher initial masses, as shown by the green and magenta 
curves in Fig. [1] In this case we typically find initial clus- 
ter masses (in stars) 15-25 times the present-day mass, for 
a Salpeter slope. Of course, the masses of proto-GC clouds 
must be even higher, depending on their star-formation ef- 
ficiency. 



4 CONTRIBUTION TO THE GALACTIC 
STELLAR HALO 

The ratio between the ejected stellar mass and the present- 
day mass and its dependence on the IMF slope is illustrated 
in the right panel of Fig. [T] The total amount of low-mass 



(< 0.8 Mq), long-lived stars (both the first and second gen- 
eration) ejected from the cluster shown by the blue lines 
corresponds to ~ (2.5 — 3.2) times the observed, present- 
day mass of globular clusters, for a Salpeter IMF in the 
case where no second generation stars are lost as assumed 
by DCM07. This amount increases for a steeper IMF. These 
stars must contribute to the population of the Galactic halo. 

Recent observations have found indications for chem- 
ically polluted, 2"^^ generation stars in the Galactic halo , 
with a frequency of fs ~ 1.4- 2.5% (jCarretta et"aLll2010bl : 
iMartell fc Grebelll2O10D . While IVesperini et al.l l|201oh have 
examined the ejection of these stars with hydrodynamic clus- 
ter models, we follow here a different approach. If these 
stars originate from the population of present day GCs, 
we can easily infer the escape fraction of 2"'' generation 
stars e^? = (1 + (1 - fv)l{fs x M^^iJ M°^))-^ , where 
Mac /Mh^io ~ 2% is t he fraction of present-day GC of th e 
total stellar halo mass (iFreeman fc Bland-Hawthornll20o3 '). 
For fs = 1 (2.5)% we obtain elf = 0.43 (0.65), i.e. a loss of 
approximately half of the second generation low-mass stars. 
With such a loss, the initial cluster masses and the total 
amount of ejected stars must be even higher than discussed 
above. The corresponding values are shown in Fig. [T] with 
green and magenta lines. Typically both initial and ejected 
mass are increased by a factor 1.7-3.5 compared to the case 
of Cll=0. For a Salpeter slope the mass of ejected low-mass 
stars is then ~ 5-10 times the present day GC mass. 

From the current total mass of halo Galactic globu- 
lar clusters of ^ 2 x 10^ Mq and the total halo 
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mass Mhaio ~ 10^ Mq jFreeman fc Bland-Hawthornll2002l ). 
i.e. tiie above 2%, we therefore find that low-mass stars 
ejected from the present-day population of GCs make up 
~ 5-8% of the mass of halo stars if e^^^O, or 10-20% (for 
the above values fs) if all halo 2"^* generation stars also 
come from these clusters. These numbers could be a fac- 
tor 0.75 lower if a lower mass-to-light ratio o f M/Lv = 1-5 
iDu bath fc Grillmair" 1997; Larson ct al. 2002, cf ) instead of 
2 adopted by (Freeman fc Bland-Hawthorn . 2002. and oth- 
ers) was more appro priate. 

For comparison. ICarretta et all (|2010bl ) estimate a min- 
imum contribution of GC stars to the st ellar halo of 2.8%, 
but up to a factor ~ 10 more, while from iMartell fc Grebell 
1120101) one obtains a contribution of ~ 17.5% 0. These 
estimates, based on the observed fraction of 2"'' genera- 
tion stars in the halo and on models accounting for mul- 
tiple stellar generations, agree with ours. Other authors, 
using calculations of the GC survival fract ion or their de- 
struction rate (cf. ICnedin fc Ostriker 1997t). estimate con- 
tributions of < 10 % jParmentier fc Gilmord 120071 ') or less 
(jBolev et al.l l2009l . 3-8 %) to the Galactic halo from the 
present-day GC population. Other estimates, based on a 
variety of different initial mass function s of st ellar clusters 
range e.g. from 4-40 % llFall fc Zhand l200ll ). ~ 40-50% 
jBaumgardt. Kroupa, fc Parmentieil 2008h . or up to 30-80 



% (|Bolev et al.l l2009l'). However, these studies neglect the 



observed multiple stellar generations in GC and their impli- 
cations. 



5 IMPLICATIONS FOR THE INITIAL 

GLOBULAR CLUSTER MASS FUNCTION 

The fact that GCs s how - now basically by definition (cf 
ICarretta et al.|[2010bl ) - stellar generations, distinct by their 
chemical abundances, allows us to take a new step in con- 
straining their initial cluster mass function. Indeed since 
stars showing the abundances characteristic of 2"^^ gener- 
ation stars have recen t ly been found in the G alactic halo 
l|Carretta et all l2010bl : IMartell fc Grebell I2OI0I ). their fre- 
quency among normal halo stars provides interesting, new 
constraints within the framework of the model examined 
here. 

Consider two limiting cases. First let us assume that 
all 2"'' generation stars found in the halo originate from 
the present-day population of GCs. As discussed above, one 
then finds for the initial mass of GCs typically Afini ~ 
(15 — 25) X Mobs- The globular cluster initial mass func- 
tion (GCIMF) must then be equal to the observed one - 
commonly described by a log- normal with a characteristic 
mass Mc ~ 1.4 x 10^ M© (cf. lHarrislll99ll ') - but with 
increased by a factor 15-20, i.e. Mc ~ (2.1 - 3.5) x lO'' Mq. 
Besides this, there is no room left for other proto-GCs, since 
these would otherwise contribute - by definition - additional 
2"^^ generation stars to the halo. 

Now assume the other extreme, namely none of the 
observed 2"^^ generation halo stars are from the observed 
GCs. This corresponds to e^=Q, and we know then that 



* Their estimated 50% of the halo mass corresponds to the total 
stellar mass including the full mass spectrum. For a Salpeter slope 
one has a fraction ^ 0.35 of low-mass stars; i.e. 17.5%. 



the present-day GC population had an initial mass func- 
tion with a characteristic mass ~ (8-10) times the present 
value (Mc ~ (1.1 - 1.4) x lO'' M©) and a total mass of 
M^i ~ (1-6 - 2) X 10* Mq. In addition, however, other, 
dissolved GCs must be invoked to explain the presence of 
these peculiar stars in the halo. Let us assume that the to- 
tal GCIMF is given by a power law with a slope /3 = — 2 , 
as often studied Ce.g. lFall fc Z hang"200ll: lBolev et al1l2009l '). 
The lowest normalisation we can chose is the one tangen- 
tial to, i.e. osculating the log-normal init ial mass fun c tion o f 
the present-day GCs, as discussed e.g. in iBolev et al.l(|2009h . 
From this we compute the total initial mass of GCs to be dis- 
solved, M°p"'^'=% by subtracting Mjjp from M*°' given by 
th e integral of t he G CIMF over the same range considered 
bv lBolev et al] (|2009l fl. Since our osculating ma ss is (8-10) 
times higher than the one of iBolev et ah! (|2009l ) we obtain 
Ml°l = (1.6 - 2) X 10® Mq, from which 90% (or more) is in 
globular clusters which must be dissolved. Since the fraction 
of low-mass (< 0.8 Mq) stars is fl% ~ 0.35 for a Salpeter 
slope, all the globulars then contribute approximately 50- 
70% of the present-day stellar halo mas^ 

From Eq. [Tjand counting the fraction (1 — /p) low-mass 
stars, we finally obtain the total amount of 2°"^ generation 
stars produced in these clusters M^"'*'"" « 0.075 x Ml°l , 
which corresponds to a fraction /s ~ 12-15% of the halo 
mass. If we adopt a lower mass-to-light ratio {M/Lv = 1.5) 
and assume that 30% of the present GC mass consists of stel- 
lar remnants, the expected fraction of 2"^^ generation stars in 
the halo may be somewhat lower, ~ 6-8 %. To compute this 
value we have implicitly assumed - in the absence of other 
information - that all GCs can be described by the same 
values of fp and as those derived from the present-day 
GCs. 

As can be seen, our theoretical prediction for the frac- 
tion fs of 2"^^ generation stars in the halo is consider- 
ably larger t han the curren t observ ationa l values of ~ 1.4- 
2.5 % from ICarretta etalj (|2010bl ) and IMartell fc Grebell 
(|2010l ). The simplest conclusion from this contradiction is 
that the GCIMF cannot be a simple power law with /3 = —2 , 
as suggested by numerous au thors (cf. FaU fc Zhanj |200ll : 
iBolev et al.ll2009l : lElmegreenlfeoiOal ). at least not over the 
mass range considered here. However, to reconcile our pre- 
diction with the observed value of fs, one would need to 
strongly reduce range of the initial cluster masses, since each 
decade in mass contains the same amount of total mass for 
this power law distribution. In other words, reducing the pre- 
dicted fs by a factor 4 or more would imply an initial cluster 
mass function over less than 1 dex, compared to our assump- 
tion of log(m) £ [3.3, 7.3]M0. Similarly, postulating e.g. that 
clusters below a certain mass (say the present day value of 
Mobs < 4. X 1O"*M0 suggested bv ICarretta et al.1 l|2010bl ')') 
will not become globulars, does not solve our problem. Alter- 
natively, in most clusters the fraction fp of unpolluted stars 
could be higher than the value observed in present-day GCs, 
in which case one could avoid "overproducing" the number 



^ Approximately over log(M) £ [3.27,7.27], corresponding to a 
range of magnitudes V from -12 to -2 for the present-day GC 
mass function. 

^ Assuming M/Ly = 1-5 instead of 2 (cf. above), this percentage 
would be lower by a factor 0.75. 
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of 2"^* generation stars in the halo. In this case, however, we 
cannot properly speak of an initial mass function for GCs, 
since these objects with much higher values of /p cannot be 
the progenitors of the present-day GC population. 

We are therefore naturally drawn to abandon the pic- 
ture of a "universal" initial power law mass function for all 
clusters, including super star clusters, young massive clus- 
ters etc. and for progenitors of present-day GCs. Then, as 
already discussed above, the observations of the 2"^* gen- 
eration halo stars can be un derstood if the GC I MF i s log- 
normal, as e.g. propos e d by IVesperini fc Zepj (|2003h and 
iParmentier fc Gilmor3 (|2005l . 120071. However, other ini- 
tial mass functions, e.g. a power-law with a turn-over or 
Schechter-type functions, cannot yet be excluded. 

In any case we have shown here that the observed frac- 
tion 2°"^ generation stars in the halo can in principle pro- 
vide very useful information on the distribution of the ini- 
tial masses of globular clusters, the GCIMF. Of course, 
our analysis does not constrain the initial mass function 
of other (non-globular) clusters. In fact, since the percent- 
age of halo stars originating from GCs is typically < 20% 
in our scenarii, there is room for other clusters, accreted 
satellites, or others to provide the rest of the present-day 
stellar halo. After many recent studies proposing a "uni- 
fied" picture for the formation and evolution of clusters of 

! Fl I I I 

all kinds inclu ding GCs (cf. Brodie & Stradcr 2006; Mcurer 
19951: 
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Fall fc Zh ang 2001 ; Vcspcr ini fc Zepf 2003; Elmegreen 



2010al . and references therein), it may well be that recent 



progress on GC stars and the finding of 2'^ generation 
stars among the halo population forces us again to revise 
this picture. One of the main questions arising now is ac- 
tually what distinguishes "normal" clusters from globulars 
and "globulars-to-become" , i.e. what causes a cluster to form 
one or two separate stellar populations. Is this e.g. related 
to their initial central density, to external conditions, or 
maybe to co mpletely different formation scenarios , as e.g. 
sugges ted bv lSearle fc ZinnI l| 19781) ; iFreemanI l| 19931 ) ; iBoketl 
(|2008l ) or can this be understood within the framework of 
curr ent hydr odynamic and cosmological formation models 
fe.g lBolev et a l. 2009; Elmegreen 2010a)? 

New detailed (hydro-)dynamic models of cluster for- 
mation and evolution taking into account recent in- 
sights gained from scenarii explaining the detailed be- 
haviour of observed ab undance patte rn in GC stars (e.g . 
Decressin. Baumgardt fc Kroupa 200^; lD'Ercole et al.l200a ; 



Decressin et al.ll2010l ') are clearly likely needed to progress 
further on this issue. In parallel it will be useful to firm 
up the first studies of second generation stars found in the 
Gala ctic halo, as they curr ently suffer e.g. from poor statis- 
tics (jCarretta et al.|[2010bl ) or from unce rtainties in obser- 
vatio nal criteria identifying these stars (|Martell fc Grebell 
l20ld ). 



4.5 




>^ Ricotti (2002) 



0.5 



1 

IMF slope 



1.5 



Figure 2. Ionising photon production of globular clusters nor- 
malised per baryon currently locked up in stars, log??', plotted as 
function of the massive stars IMF slope x. Blue, green, and ma- 
genta lines show respectively values of r)' for IMF slopes allowing 
one to reproduce the observed fraction /p (boundaries shown by 
dotted lines) as suming e£^ = 0, 0.43, and 0.65 respectively. The 
value of r) from iRicottil l|2002l ') computed for metallicity Z = 0.03 
Zq and the Salpeter slope (x = 1.35) is shown for comparison. 



put (i.e. the total number of photons emitted above 13.6 
eV) during the life of a GC normalised to its current num- 
ber of baryonfl, 77', as a function of the IMF slope for the 
dynamical scenario described previously. The Lyman con- 
tinuum flux was compu t ed using the evolutionary synthesis 
code of lSchaereil (|2003l l; iRaiter. Schaerer fc Fosburvl (|2010f l 
for a low metalhcity Z = 1/20^© typical of GCiQ. 

Interestingly, the predicted ionising photon output is 
quite independent of the IMF slope (see Fig. [21), since both 
the metals explaining the observed abundance pattern and 
the ionising photons are made primarily in massive stars. 
Should the "pollution" of 2"'^ generation stars be related 
to massive AGB stars, we expect a similarly large emis- 
sion of ionising photons per baryon, since the ratio between 
the initial and present-day mass is similar in both scenar- 
ios (cf. above). In the AGB scenario the dependence of ry' 
on the IMF should, however, but somewhat stronger, since 
the masses responsible for the stellar ejecta and the ionising 
photons are more distinct than in the fast rotating, massive 
star scenario. Of course, one also finds an increased output of 
ionising photons per unit present-day mass (or baryon num- 



6 THE CONTRIBUTION OF GLOBULAR 
CLUSTERS TO REIONISATION 

Since the initial masses of GCs may be substantially larger 
than their present day values, their output of ionising radia- 
tion and their contribution to cosmic reionisation also needs 
to be revised. 

In Fig. [2] we plot the predicted H ionising photon out- 



To convert photon/baryon into photon/mass one has e.g. 
IMq = 1.19 X 1057 baryon. 

* Adopting a different metallicity leads to small changes (typ- 
ically < 0.1—0.2 dex). Similarly, using the stellar evolutionary 
tracks of fast rotating stars used in the study of DCM07 would 
lead to relatively small changes compared to the effects discussed 
here. 
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ber) if loss of second generation stars is allowed (green and 
magenta lines) compared to the case where Bll ~ (blue 
lines). This increase is simply due to higher initial mass of 
GC s, discus s ed pre viously. 

iRicottil (|2002h has estimated the output of GCs per 
baryon as r/' = r; x fdi, with 77 ^ 900(0 and fdi ~ 2-10 
(with a maximum of fdi < 100), where fdi - the equivalent 
of our ratio Mini/Mobs - is a factor accounting for dynami- 
cal disruption of GCs during their lifetime. In our case, the 
total photon output is log(r;') ~ 4.8-4.9 photon/baryon if 
all second generation stars remain within the cluster, which 
is approxi mately a factor 5 higher than the typical value 
adopted bv lRicotTH (|2002l ') for fdi ~ 2. The amount of emit- 
ted ionising photons must be even higher if 2"^* generation 
stars were lost from these clusters (gll > 0), as illustrated 
by the green and magenta lines. 

Using simple but elegant arguments to estimate the 
number of ionising photons emitted per baryon during the 
formation of GCs, and using a Press-Schechter model to 
compute the GC formation rate, IRicottil l|2002h has esti- 
mated the contribution of GCs to cosmic reionisation. He 
has shown that the old GCs produced enough ionising pho- 
tons to reionise the inter-galactic medium at 2: ~ 6, if the 
escape fraction of ionising photons, /esc, from these objects 
was close to unity. Our finding of a high ionisi ng phot o n pro - 
duction in GCs reinforces the conclusions of IRicottil (|2002l ) 
and leaves room for lower values of /esc < 1 or for other 
less favourable assumptions (e.g. uncertainties in the age of 
GCs). It appears that old GCs, formed as massive super star 
clusters shortly after the Big Bang, could provide a signifi- 
cant if not dominant source of UV radiation to reionise the 
Universe at high redshift. In any case, it seems unavoidable 
to seriously consider their contribution. 

With a typical mass of ~ 1.4 x 10^ M0 a young proto- 
GC is expected to have a peak UV (say 1500 A) magnitude 
of Mab —16.7 at an age '--^2 — 3 Myr, before fading 
by ~ 4 mags within ~ 20 Myr just due to stellar evolu - 
tion (see e.g. models of lLeitherer et al1ll999l : ISchaereij|2003l ) 
or faster when evaporation (mass loss) sets in. At redshifts 
of z ~ 7-10 this would correspond to a typical UV rest- 
frame magnitude of mAB ~ 30.3 at the peak, just slightly 
fainter than the current detection limits of the deepest near- 
IR images taken wit h the WFC3 camera onboard HST (cf. 
iBouwens et al.|[20ld ). In any case, single, massive proto-GCs 
during their youth might be detectable in situ with current 
instrumentation, and are certainly well within the reach of 
even deeper near-IR observations, which will be achieved 
with the JWST. However, whether the two proposed sce- 
narii (massive AGB or fast rotating, massive stars as the 
origin of the bulk of material out which 2""^ generation stars 
form in GCs) can be distinguished observationally at high 
redshift, appears a priori quite difficult, if not impossible. 



7 SUMMARY AND CONCLUSIONS 

In light of the recently recognised, general existence of mul- 
tiple stellar generations in globular clusters (GCs) implying 
significant losses of 1°' generation stars from these clusters, 

^ Making the same assumptions as Ricotti we confirm this value. 



we have re-examined the initial masses of GCs, the contri- 
bution of low-mass stars ejected from GCs to the stellar halo 
of our Galaxy, and the contribution of GCs to the ionising 
photon production necessary to reionise the inter-galactic 
medium at high redshift. 

These quantities have been estimated from the chemical 
and jiynamical model of iDecressin. Charbonnel &: Mevnej 
(2007), which successfully reproduces the main obser- 
vational constraints from 1*^* and 2"^* generation stars, 
by invoking pollution from fast rotating massive stars. 
The main free parameters of this model are the slope 
of IMF for high masses (> 0.8 M©, the IMF being 
fixed to the observed log-normal distribution for lower 
masses), the relative number of l'''/2"'* generations stars, 
given by the fraction fp — 0.331!q q8 generation 
stars deter mined from t h e deta iled spectroscopic obser- 
vations of C arretta et al.l (j2010bl ). and a dilution param- 
eter d ~ 1.15 i nferred from the Li-Na anticorrelation 
observed in GCs jPecressin. Charbonnel fc Mevnetl l2007l : 
ICharbonnel fc DecressinI 201ll ). 

The dynamical scenario we have explored allows for the 
evaporation of stars from the 1^^ generation (correspond- 
ing to an escape fraction of 2""^ generation stars of zero, 
^LL =0), or from both generations, as suggested by recent 
observations finding stars characteri stic of the 2°*^ genera- 
tion in GCs in the Mil ky Way halo |Carretta et al.ll2010bl : 
iMartell fc Grebelllioioi ') . The latter case translates to elf ~ 
0.43-0.65. 

We have obtained the following main results for an IMF 
with a Salpeter slope above 0.8 M0: 

• The initial stellar masses of GCs must have been ~ 8-10 
times larger than the current (observed) mass, when no sec- 
ond generation stars are lost, in agreement wit h the earlier 
results ofjPrantzos & Charbonnel (2006) and Ca rretta et al.l 
(|2010bl ). If all 2"^^ generation halo stars originate from the 
present population of GCs, the initial cluster masses must 
have been ~ 25 times larger than the current mass. 

• The mass in low-mass stars ejected from GCs must be 
~ 2.5 — 3.2 times their observed, stellar mass if all 2"^^ gen- 
eration stars were retained, or ~ 5-10 times the present 
day mass if ~ 0.43-0.65. These numbers translate to a 
contribution of 5-8% or 10-20% respectively of the ejected 
low-mass stars to the Galactic stellar halo mass. We have 
compared our estimate with earlier values obtained from 
various methods (cf. Q. 

• The observations of 2"^^ generation stars in the Galac- 
tic halo can constrain the initial mass function of the GC 
population (GCIMF). In particular we have shown that a 
power-law with a slope ~ —2, as often assumed, is in con- 
tradiction with recent determinations of the fraction of 2"^* 
generation stars in the halo, whereas a log-normal GCIMF 
is compatible with these observations. This finding revives 
the question about a common mass function and about the 
physical processes leading to a distinction between globular 
clusters with multiple stellar populations and other clusters. 

• Due to their high initial masses, the amount of Lyman 
continuum photons emitted by GCs during their youth must 
have been substantial. Indeed, we find that their output cor- 
responds to a total number of ionising photons emitted per 
baryon, tj' ^ lO^-S-^-^ for elt = 0, or ~ 1.7-3.5 times more if 
6ll ^ 0.4-0.6. Our results reinforce the conclusion of lRicottil 
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that GCs should contribute significantly to reionise 
the IGM at very high redshift (2 > 6). Individual, young 
proto-GCs with typical masses few times ~ l(f Mq could 
just be detectable at high redshift in ultra-deep images with 
the HST, and are certainly within the reach of the JWST. 

The dependence of the initial and ejected masses on the IMF 
slope has been illustrated in Fig. [T] The ionising photon 
production is found to be quite insensitive to the high mass 
IMF, since both the ejecta "polluting" the 2"^* generation 
stars and the Lyman continuum flux originates from massive 
stars. Our main results should also be valid for the massive 
AGB scenario, at least qualitatively. 
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